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ABSTRACT 

A system was designed to investigate experimentally the mechanism 
of nucleate pool boiling heat transfer using commercial grade liquid 
nitrogen as a working fluid. aA circular horizontal flat plate five square 
centimeters in area was utilized as a boiler surface, 

The effect of various surface parameters on the characteristic heat 
flux versus (T,-T,) curve at atmospheric conditions was determined. 
Surfaces used included a highly polished mirror surface, a milrror surface 
coated with oxide, a mirror surface coated with grease, a roughened 
surface, and a roughened surface coated with Teflon, all fabricated from 
commercial electrical tough pitch cope. The final surface used was a 
Nickel 200 highly polished mirror surface, 

The data from the copper surface with a mirror finish was in agreement 
with the data of previous investigations. The effects of roughness, 
contaminants, hysteresis, and boiler surface material were discussed. A 
comparison was made on incipient boiling heat fluxes, The results indicate 
that the surface conditions play a major role in nucleate pool boiling 


heat transfer to liquid nitrogen, 
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SECTION I 
INTRODUCTION 
= ohgrowed Background 

Recent technological advances in the space program have brought 
forth great improvements in the ability to obtain and use cryogenic fluids 
for modern research and industry. Because of its infancy, however, there 
is still much left unexplored in the realm of cryogenic heat transfer, 
The uses for the liquified gasses are as varied as their properties. For 
instance, liquid nitrogen, because of its relative ease in handling and 
chemical properties, has been used in quick-freeze food preparation, 
while research laboratories use it as a combination refrigerant-heat 
sink, The space program has made extensive use of this coolant as an 
inert pressurizer. Photographic techniques have been refined using liquid 
nitrogen, Even the medical profession has found uses for liquid nitrogen 
in blood and tissue preservation, and in cryosurgery. 

All of the above uses result in heat being transferred to the liquid 
nitrogen, As observed in Figure 1, the heat flux (Q/A) and temperature 
difference (AT) determine the mechanism of heat transfer, whether it be 
by convection, nucleate or film boiling. Because most uses involve the 
nucleate boiling regime, it is this region that is the general subject 
of this investigation. 

oo PE loud ast ire. Previous Research 

Previous investigations dealing with nucleate boiling using non- 
cryogenic fluids have attempted to relate effects of surface conditions, 
among which are included roughness, effects of contaminants, material, 
and material history. Historically, Corty and Foust [8] were among the 


earliest to attempt correlation of the effect of surface roughness on 
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nucleate boiling heat transfer coefficient, using diethyl ether and n- 
pentane. . Generally, their work showed that the polishing of a copper 
and nickel test surface affected both the position and slope of the 
boiling curve. In 1962, Berenson | 4 | furthered the state of the art by 
comparing nickel, copper and inconel, with n-pentane using more refined 
techniques. Others such as Westwater, Clark and Strenge [23] were 
experimenting with boiling water and high speed*motion picture studies 

of active sites from scratches and pits in order to determine the nucleate 
boiling mechanism. Surface contaminants were the subject of investigations 
by Averin (oil coating) { 1] » while Young Satan) [26] » uSing water, 
presented data on the effects of Teflon coatings on a roughened surface. 

A theoretical investigation was carried out by Bankoff | 3 | » from 
which he concluded that nice oe from a homogenous liquid was nearly 
impossible, but he postulated that nucleation proceeds from a pre-existing 
vapor phase by analyzing fluids in wetted and unwetted cavities. This 
postulate is verified by Griffith and Wallis [10] in their classic 
experiment. Bankoff concluded that for wetted cavities, vapor embryos 
will collapse completely and the cavities will be nullified as nucleation 
centers. He further states that if the surface is only partially wetted, 
there will be some steep walled cavities which retain vapor phase, A 
slight increase in temperature or decrease in pressure would cause the 
nucleus to grow spontaneously. 

Even with all these efforts, the exact mechanism of boiling and the 
effects of various parameters are Fgeenore ere understood. A number of © 
investigators have attempted to correlate data. One of the equations 
which better accounts for surface variables is that of Rohsenow [19] . 

In his equation, C.¢ is an empirically determined constant for the fluid 


surface combinations, indicating that for a theoretical evaluation there 
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is a need to point out variables with each and every surface combination, 
For example, amounts of dissolved and/or adsorbed gasses, impurities, 
contaminants, and the nature of the heating surfaces may vary from one 
investigation to another, making correlation extremely difficult. 

There has been very little reproducible data presented considering 
the effects of surface conditions on boiling heat transfer using cryo- 
genic fluids. Summaries of data presented by Seader and Miller { 21] , 
and, more recently, by Brentari, Giarratano and Smith | 7 | » indicate a 
very wide spread in data due to surface conditions alone. It is difficult 
to correlate experimental results with theory and the variables are hard 
to separate. This has resulted in many investigations which have been 
undertaken without control of individual variables. 

Ob jective 

The objectives of this research were fourfold: first, to design a 
suitable test device for determining heat rate (Q/A) versus temperature 
difference (T,-T,) for liquid nitrogen and liquid helium, using as much 
existing equipment as was consistent with the objectives and budget; 
second, to establish experimental procedures which would yield reproducible 
nucleate pool boiling data for a given set of boiler surface conditions 
with liquid nitrogen; third, to make provision for photographic analysis 
of bubble nucleation using a high speed camera; fourth, to take data on 
five copper surfaces to include a highly polished mirror finish, a mirror 
finish coated with oxide, a mirror finish coated with grease, a roughened 
surface, and a roughened surface coated with Teflon. 

In order to design the equipment to accomplish the objective, the 
following conditions were established: 

1. The working fluid was liquid nitrogen because of its availability 


in high purities at low cost. In addition, the physical properties were 
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well tabulated and liquid nitrogen exhibited relative ease in handling. 

2. The system was open cycle in that nitrogen was allowed to boil 
aways; however, contamination of the boiler and working fluid from 
condensing oxygen and moisture was to be minimized, 

3. The boiler surface was a horizontal flat plate, because it 
offered the simplest geometry to analyze, and data was available for 
comparison, Additionally, this design permitted ease in the regulation 
of test surface finishes and in changing test surfaces. 

4, Provision’ was made for visual observation of the boiling surface 
during data collection so that future studies could incorporate high speed 
motion picture analysis. 

5. The system was compatible with liquid helium as a working fluid 
for subsequent investigations. 

6. The system was as safe, simple, and flexible as possible. 

With these conditions uppermost in mind, the final design resulted in the 
apparatus shown schematically in Figure 2. A photograph of the actual 
system is given in Figure 3, Figures 4, 5, and 6 show an exploded view 
of the boiler-heater assembly, a sectional drawing of the boiler test 
sections, and an assembled boiler test section. Figures 7 and 8 are a 
schematic drawing of the vacuum control system and the temperature sensing 
system, General cryogenic design background was provided by Scott [22] 


and White [25] ° 
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SECTION 2 
DESCRIPTION OF EQUIPMENT 
pi . Boiler Test Section 

(1) Heater, . An electric heater (Figure 4) was fashioned from a 
19 inch long inconel-sheathed (.045 inch diameter) Nichrome wire heating 
element, manufactured by the American Standard Company. Insulation 
material between the sheath and the wire was magnesium oxide, This heat- 
ing wire was wound by hand around a 1/16 inch starting mandrel,. Final 
diameter of the button~shaped heater was slightly less than one inch. 
Resistence of the heater element was 42 ohms. 

Previous experimenters studying boiling heat transfer have had a 
choice between control of heat flux or control of temperature difference 
in their investigations; similarly, the choice occurred here, While the 
control of temperature difference allows a wider range of operation over 
the characteristic boiling curve, with ease in operating in the transition 
region (Figure 1), the increased complexity of the apparatus for this 
initial study, plus the fact that only the nucleate region was under study, 
justified the use of an electric heater to control heat flux. Data of 
previous investigations [ 7 | indicated that for a variety of surface 
materials, expected heat fluxes on the order of a maximum of 20 watts/om” 
were needed in order to attain the maximum nucleate boiling flux, A 
boiler surface ofl fine cm* was selected as being compatible with this flux 
stipulation, and available space and heaters. In order to maintain one- 
dimensional transfer, the heater had to have a dimension very close to the 
maximum diameter of the boiler. 

For the expected heat flux, maximum voltage across the heater was 
calculated to be 65 volts, with a maximum current of 1.55 amps. The watt 
2 


density was calculated to be 37.8 watts/inch” per inch of heater, conform- 
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ing to standard practice. A density of 45 watts/inch® per inch of heater 
was considered as a maximum for this type heater, A sheathed heater was 
utilized to reduce any electrical noise in the vicinity of the test surface 
and to provide a surface capable of withstanding the silver solder braze. 
(8) Transition Section. The cylindrically shaped transition 

section of diameter .9935 inches and length .315 inches (Figure 4) was 
placed between the heater and the nickel boiler or between the heater and 
disc in the case of the copper boiler. The function of this section was 
to even out any heat flux discontinuities from the heater, which was 
brazed to it using silver solder, 

wey sSinless Seoel iso, < The Stainless Steel 321 disc shown 
in Figure 4, had a thickness of .025 inches with a diameter of 1.00 
inches, An exact match to the diameter of the transition section and 
boiler was compromised, as the disc used could be punched out of sheet 
stock rather than turned down from larger stock. This disc, used only 
with the copper boiler surface, provided a thermal resistance, insuring 
that operation of the heater attached to the copper sample was at the 
Same temperature as tien it was attached to the Nickel 200 boiler, 
This disc was placed between the transition section and boiler test 
surface and soldered in place. 

(4) Boiler Plate. The boiler plate (Figure 4) was cylindrical in 
shape and was made according to the dimensions shown in Figure 5, The 
length of each cylinder was selected to give a temperature distribution 
large enough to determine heat flux density with an error of less than 
ten percent at maximum heat flux. (See Appendix C for sample calculation.) 

(5).Boiler Diaphragm and :Enclosure. The boiler diaphragm, used to 
support the boiler on the enclosure, was .020 inches thick and diameter 


of 1.78 inches, The enclosure was a hollow capped cylinder of 1.78 inches 


16 


O.D. and wall thickness of .035 inches. Both items were made from stain- 
less steel and are pictured in Figure 6. The diaphragm was 
manufactured with a cylindrical lip .035 inches thick an: ,125 inches 
long which fit snugly inside the enclosure and was used to secure a leak- 
tight braze with relative ease, 

The two factors determining the thickness of the diaphragm were 
strength and high resistance to thermal losses. A compromised dimension 
of .020 inches provided safety and good thermal characteristics. The 
heater enclosure was designed for low stresses and to allow room for the 
associated instrumentation and power wires, 

SY mene Dysarts Vacuum Systems 

For this investigation two vacuum systems were needed, A medium 
system (107mm Hg) was utilized for thermal insulation of the test boiler 
and inner dewar, while a low vacuum system (1073nm Hg) was used for 
environmental control of the boiler fluid, Both are shown schematically 
eT aeirer2 | 

(1) .Thermal Insulation System, A large capacity Welch Duo-Seal 
1402-B Two-stage Mechanical Pump served as a roughing pump for the system 
and fore-pump for the diffusion pump. Rated capacity was 140 liters/min- 
ute of free air displacement with an ultimate vacuum of 1074mm Hg. The 
pump was vibration-isolated from the pumping unit, using a mounting 
platform on four small jacks extending to the floor. A Consolidated 
Vacuum Corporation MCF 300 Fractionating type Diffusion Pump provided the 
medium vacuum, The rated capacity was 290 liters/second at 0.4 micron 
with an ultimated vacuum of 2. x 107’mm Hg. Dow-Corning 705 Silicon Oil 


was used in the diffusion pump because of its resistance to oxidation 


and very low backstreaming characteristics. 
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Safety features installed for diffusion pump protection included a 
thermal cut-out switch, mounted on the cooling coils of the diffusion 
pump, and a diaphragm vacuum microswitch, located in the fore—pump line. 
A liquid nitrogen cold trap was placed between the diffusion pump and 
the test unit to capture condensable vapors and reduce backstreaming, A 
15 inch 0.D. x .035 inches wall stainless steel pipe connected aheid: ri. 
trap to the valve manifold, From the valve manifold, the pipe diameter 
decreased to 4 inch 0.D. An Anaconda copper and Mason-Renshaw stainless 
steel flexible ~ connector were attached to quick-disconnect couplings 
leading respectively to the dewar and the test unit to provide ease and 
speed in assembly. 

Wherever possible, stainless steel was used to reduce the problem of 
outgassing within the system. Standard copper tubing elbows were used in 
conjunction with the stainless pipes, Viton ‘Af? O-rings were used at all 
flanges because of their low outgassing characteristic. Thermocouples 
and power leads for the boiler unit penetrated the vacuum system through 
a Conax MHC Multiple-hole Saoemne Gland (ceramic followers with a Teflon 
sealant). This packing gland assembly and the copper elbow fittings were 
brazed to the stainless pipes using silver solder. 

The electrical design for this vacuum system control is show in 
Figure 7, The safety features prevent damage from loss of water, vacuum 
or electrical power, Loss of water pressure, dangerous in that pump oil 
could reach the flash point, would cause the diffusion pump wall tempera- 
ture to rise to above 250°F and thereby open the hold-down control coil 
circuit and shut down the pump. Similarly, a loss of fore—pump vacuun, 
which would cause oxidation of the diffusion pump oil, would also open 


the circuit. This feature also prevented energizing the pump without 
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sufficient vacuum, A loss of electrical power would also open the hold- 
down control coil. In all cases of automatic shutdown, manual restart 
would be required. 

The thermal insulation vacuum pumping system with control power was 
assembled as a separate pumping unit for ease in control, access, and 
flexibility (Figure 3). With minor modifications, the system could ke 
used as a portable system in other studies, A vacuum on the order of 
8 x 107Smm Hg or better, attained with relative ease, was used for test 
runs. 

(2.)' Environmental: Control. The mechanical pump for this system, 

a National Research Corporation type 2S, although rated at only 33 liters/ 
minute, proved very satisfactory, The ultimate vacuum of this pump was 

3 x 1074mm Hg. A liquid nitrogen cold trap was provided between the 

pump and the dewar. This vacuum system provided control of the environ- 
ment in the inner dewar. Prior to a run this pump was used to outgas 

the test surface and dewar. In future studies, it can be used in studying 
the effect of reduced pressure on nucleate pool boiling. 

, ees Dewars 

The inner dewar (I.D. = 3.1 inches), shown schematically in Figure 
2, held the working fluid which surrounded the boiler assembly, and was 
provided with an O-ring and flange for vacuum seal in order to control 
pressure over the fluid. This dewar was connected to the insulation 
vacuum system by means of a 4 inch Pyrex-Kovar Housekeeper Scal which was, 
in turn, connected to the quick-disconnect coupling. A dewar with a 
permanent" vacuum was not used, as helium penetrates warm pyrex glass at 
a slow rate, thereby losing the insulating effect. 


The outer dewar (I.D. = 4.7 inches), filled with liquid nitrogen and 
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used as a radiation and intermediate thermal barrier, was manufactured 
with a sealed permanent vacuum of 107%mm Hg. The dewars were each 
fabricated with a double silvered mirror radiation shield and an off-set, 
aligned, unsilvered, vertical strip 0.8 inches wide, which allowed for 
visual observation of the boiler surface during boiling. 

This double pyrex dewar arrangement provided the needed containers 
for controlling the environmental conditions affecting the working fluid. 
At no time was boiling ever observed from the heater assembly due to 
radiation energy penetrating either through the unsilvered strip or 
plastic viewport in the dewar flange as long as both dewars were filled. 

ae, Instrumentation 
(L) Temperature Sensing, ISA type T copper-constantan thermo- 
couples sheathed with Stainless Steel 304. (sheath 0.D. = .041 inches, 
wire diameter = .007 inches), were selected with special limits of error. 
These high purity sensors were fabricated five feet long with a grounded 
junction, An x-ray view of this junction is pictured in Figure 9. Care- 
ful attention was given to the selection of good detectors. The tempera- 
ture difference across the copper boiler test section was expected to be 
a maximum of 8 K. Accuracy of 0.2°K was desired and space was restricted. 
Even though iron-constantan thermocouples had a larger electromotive 
force per degree, a greater degree of accuracy could be attained (Dike 
{2]) with the copper-constantan combination using the calibration 
methods outlined in Appendix A. 

After calibration, four thermocouples were distributed over the 
length of the test element (Figure 5). Each was coated with Apiezon "Nt" 
vacuum grease, then placed in a .042 inch diameter hole which had pre- 
viously been filled with this grease. Immersion depth was .37 inches. 


For this immersion depth, a conduction error of less than .04°K was 


oO 
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expected | 2]. The grease provided good thermal contact {18} with the 
boiler and facilitated the element assembly. 

The thermocouple leads were fed through the vacuum support tube and 
the Conax gland, then attached directly to a rotary switch, eliminating 
associated small errors of extension leads. A liquid nitrogen reference 
junction was used to minimize error, and the fluid was agitated in order 
to prevent stratification. A fifth thermocouple was located in the boiling 
liquid nitrogen to detect changes in bulk temperature during operation. 

The rotary switch and thermocouple junctions were located in a zone 
box (fabricated from 5 inch plywood), which prevented any undesirable 
Seebeck emfs. The thermocouple output was amplified by a Hewlett Packard 
Differential Amplifier, Model 8875 A, and fed to a Hewlett Packard auto- 
mated D.C. Digital Voltmeter, Model 405C. Amplification was necessary 
in order to obtain three significant figures from the microvolt thermo- 
couple output, as the minimum voltmeter sensitivity was 0.001 volts. 

Care was taken to insure that leads were electronically shielded 
from external noises to prevent any stray current error (following guidce- 
lines outlined by Morrison [15] ). An internal signal with variable rate 
control in the voltmeter keyed a Hewlett Packard Digital Recorder, Model 
516 B, which printed the visually displayed digital reading on three 
inch strip paper (data sheet), A schematic diagram is shown in Figure 
7. For calibration of the amplifier-digital voltmeter combination see 
Appendix A. 

(2) Pressure. Sensing. Pressure in the vacuum system was monitored 
at low vacuums with National Research Corporation, Model 724, matched 
thermocouple gauges and at high vacuums with a NRC, type 724 Cold Cathode 
Ionization Gauge Control with sensing head 524, While taking data ona 


particular surface, only the cold cathode gauge was monitored. 
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SECTION 3 
EXPERIMENTAL PROCEDURES 
Preparation of Test Surfaces 

Because of the complex dependence on surface properties, surface 
preparation, past history of the boiler, etc., the number of active 
sites of a given size cannot be predicted either theoretically or experi- 
mentally. It is hoped, however, that by describing the surface prepa- 
ration, the microstructure is sufficiently well defined for use in future 
duplication or correlation. 

A machined boiler test surface (see Figure 4) was brazed to a 
stainless steel diaphragm while, at the same time, a heater was brazed 
to the copper distribution plug, both with silver solder (m. p. 1170°F) 
and ruby flux. Upon cooling, the two partial assemblies were bonded 
together with Eutecrod 157 solder (m. p. 425° F) (95% tin, 5% silver) 
with the stainless steel disc used between the copper boiler surface and 
the distribution plug. The purpose of the lower temperature bond was to 
allow for separation of the heater for use with another test boiler 
surface, The surface of this entire assembly was then cleaned to remove 
flux, checked for continuity in the heater, and polished. A complete 
assembly is shown in test position in Figure 6. 

(1) -Mirror: Finish. The boiler surface, protruding 0.015 inches 
above the stainless steel diaphragm, was sanded by pressing against a 
rotating dry 320 grit carborundum belt. The sample was then rotated 90° 
stroked by hand in one direction only over 0 emery paper, rotated again 
90°, and stroked over 3/0 emery paper, both emerys being dry. Only as 
many strokes were used as were necessary to remove the marks from the 


previous process. Wet polishing on three 8-inch Buehler metallurgical 


22, 


polishing wheels followed. Wheel one was covered with canvas. Rotating 
the sample 90° from the previous sanding direction, and using 500 grit 
carborundum, all previous marks were removed. Visually, the surface had 
the appearance of an old mirror with many scratches showing. The second 
wheel, covered with felt-like material (Kitten Ear), was used with a one 
micron alumina heterogencous solution, The final wheel was covered with 
a softer velvet-like material (Reyvel). The abrasive used with this 
wheel was 4 micron ganma alumina and produced a fine mirror finish with 
no visible scratches. After polishing with each wheel, the test surface 
was washed using a tap water jet. In addition, upon completion of 
polishing with the third wheel, the surface was flushed with methanol 

to prevent spotting, then rotated under a hot air jet, drying the surface 
completely, Scotch Protective Tape, Type 343, with a low adhesive 
coefficient was placed on the surface to protect it from scratching and 
oxidation during installation, 

(2) Grease Coated Finish, .. A clean, mirror finished surface (as 
outlined above) was coated with Apiezon "N' vacuum grease by applying a 
small amount over the entire surface with a cotton swab. Small ridges 
were visible, however, Application of a hot air jet melted the grease 
to a thin, evenly distributed layer of estimated .002 inches thickness, 

A clear, undistorted reflection could be seen in the surface upon 
cooling. 

(3). Oxidized: Finish,-' <A clean, mirror finished surface was oxidized 
by exposing it to heat. aA propane torch, at about 200 - 230°F, was used 
to apply the heat until an opaque film covered the surface. Despite this 
treatment, the test surface's appearance was still smooth, and a reflection 


could be observed, 
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(4) Roughened Finish, A clean, mirror finished surface was again 
used as the basic surface to modify. Light taps from a small peening 
hammer were used on 5 grit flint paper which covered the surface, These 
nucleation cavities were evenly distributed over the entire surface, 

The surface was then swabbed with reagent-grade methanol and dried. The 
maximum cavity size was approximately .003 inches diameter, 

(5) Teflon—CGoated Roughened Finish, After having boiled off the 
roughened surface, it was cleaned and one coat of green DuPont 851-204 
Teflon finish was brushed on. This was wiped into the cavities using a 
soft paper towel. A second coat was immediately applied and wiped off, 
apparently filling all nucleating sites, The Teflon was then cured at 
a temperature of approximately 390 to 420°F by using a hot air jet for 
ten minutes. A slight change in hue indicated that the curing temperature 
had been attained. 

Upon cooling, the surface was polished for one minute, smoothing 
out any major asperities in the Teflon. The surface was tested for and 
found to be non-wetting for water, A practice surface treated similarly, 
however, was found to be wetted by liquid nitrogen, The maximum diameter 
of the filled cavities was approximately .003 inches, determined by the 
use of a metallurgical specimen microscope, Areas between the Teflon- 
filled cavities were highly polished, 

Testing Procedure 

Each thermocouple was inserted into its respective hole with the 
Apiezon "N't grease, wrapped a half-turn around the boiler plate at the 
same distance from the boiler surface, then secured with Teflon strips. 
The heater power wires were attached, insulated, and fastened, using the 
wire itself, to the transition section to prevent breakage. Following, 


the wired boiler test section was pushed without difficulty, into the 
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enclosure with a twisting motion to coil the wires. The upper lip of the 
diaphragm was then soldered to the enclosure with Eutecrod 157 (Figure 6). 

Upon completion of the above assembly of the boiler test plate and 
enclosure, the system was electrically re-checked for continuity, The 
boiler insulating vacuum system was connected and pump-down commenced. 
The boiler test surface was repolished just prior to a test run witha 
portable, number three wheel using gamma alumina to remove light tape 
marks, recleaned with tap water, flushed with methanol, then dried under 
a hot air jet. aA small amount of Apiezon "N" grease was placed around the 
test surface-stainless steel interface with a hypodermic needle to reduce 
"roguc't nucleation sites (Figure 6). The test dewar was then hung in 
place, connected through the quick-disconnect seal to the vacuum system 
and evacuation began, A plexiglass view port was placed on its O-ring 
located in the flange, and evacuation began inside the dewar. This 
provided an easy way to rid the dewar of moisture from a previous run 
and, at the same time, outgas the system, Thermal system pump-down 
time to 2 x 107°mm Hg was approximately one hour, During this time, the 
thermal radiation barrier dewar was cleaned and put in place with 
unsilvered strips aligned. 

Upon reaching 2 x 105mm Hg, the test dewar was vented and a small 
quantity of liquid nitrogen was admitted to cool down the test specimen 
and dewar wall. Afterward, a small amount was admitted to the outer 
dewar. Subsequently, both dewars were filled slowly in order to prevent 
thermal shocks. 

The heater was then energized to a heat flux of 5-6 watts/cm* and 
boiling commenced for ten minutes in order to degas the liquid nitrogen 
charge. The boiler was then allowed to cool. When most of the cavities 


had ceased nucleating, a quick vacuum, using the environmental pump, was 
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pulled momentarily ta.oause~ violent boiling by lowering Ts. By resuming 
atmospheric pressure conditions, the pressure. increased: thereby collapsing 
the remaining active nucleation sites, A waiting period followed, during 
which thermal equilibrium was established, This was determined by 
monitoring the boiler and bulk fluid thermocouples for steady readings. 

The barometric pressure and ambient temperature were read prior to 
a run for saturation temperature correction, 

A data collection run consisted of applying power at approximately 
10 watt intervals up to a maximum of 82.5 watts, which corresponded to 
a heat flux of 16.5 watts/cm*, and then decreasing the power at different 
intervals. Change in nitrogen height was recorded in order to determine 
the change in saturation temperature, Temperature emfs were recorded for 
the four thermocouples after changing heater power and waiting for thermal 
equilibrium to be re-established, In all cases, the thermocouple closest 
to the surface was re-read in order to duplicate the original reading 
within LW before changing the flux level, Approximate waiting time for 
equilibrium after a power change was two minutes. Thermocouple response 
to a step change was on the order of one second, while the response of 
the rest of the temperature sensing system was considerably shorter, 
depending on scaling time by the digital voltmeter, Usual time for taking 
thermocouple readings was thirty seconds at any specified heat flux. 

The system was shut down by de-energizing the heater, removing the 
outer dewar, and disconnecting the thermal vacuum, The test dewar was 
removed, allowing the boiler assembly to then warm to room temperature 
by natural convection. The thermal insulation vacuum was maintained to 
the boiler/boiler enclosure during warm-up to remove any cryopumped 
vapors in the boiler enclosure and support tube. During operation, this 


cryopumping effect helped the vacuum system such that by the end of a run, 
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system vacuum was 6 x 10° 
a Treatment of Data 
At any level of heat flux the emf values in microvolts were zero 
corrected, then converted to temperature using the working curves gene- 
rated for each thermocouple. These temperatures were plotted versus 
boiler plate axial thermocouple position to obtain a temperature distri- 
bution, A straight line always resulted. Using the slope of this line 
and an average thermal conductivity (Figures10 and 11), heat flux could 
be calculated directly, using the Fouricr heat conduction cquation: 
Q/A =k AT/&X, The wall temperature was determined by extrapolating the 
straight line tempcrature distribution to the boiler surface, 
Saturation temperature was determined from the corrected barometric 


pressure and working fluid height, using the Clausius Clapeyron equation, 


A sample calculation is given in Appendix C, 


Generally, data was processed in the metric system, as most cryogenic 


“ 


investigations make use of this system. 
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SECTION 4 
RESULTS AND DISCUSSION 

Six surfaces were evaluated. Heat fluxes as high as 16,5 watts/cm* 
were attained without complications. Some difficulty was encountered 
from ttrogue" sites between the boiler interface and stainless diaphragm 
as they tended to obscure the boiler surface, preventing photographic 
records, but the mechanism of active sites spreading on the test surface 
was observed and estimates of nucleation coverage could be made, The 
Apiexon "Nt grease applied, after a preliminary test, at the interface 
reduced the problem eonstaeiuy: The effect is considered minor, as the 
temperature distribution was linear over the four thermocouples. Any 
heat flux distortion would have been reflected in the reading of the 
thermocouple closest to the surface, 

A comparative summary of each surface is given in Table I where 
AQT equals the difference between the boiler surface temperature and the 


Saturation temperature at liquid pressure. 


TABLE I - NUCLEATE BOILING DATA SUMMARY : 


Surface Incipient | aT Full Surface| AT Last Observed AT 
Condition | Boiling Nucleation Nucleating Site 
Flux w/em Flux w/em on pecrce ang 
Flux w/cm2 


Mirror 4 

Copper #2 

Nickel 

Oxidized 


we far 
es | Loy 
me | a2” 1s) ae 


* Indicates very smal] heat flux and corresponding AT. 
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Typical data points with estimated error limits for an electrical 
tough pitch copper mirror surface are given in Figure 12. A curve showing 
the experimental results of Roubeau [29] , who used a copper boiler of 
similar size and oriented surface, is shown for comparison. The 
theoretical correlations of Rohsenow and Forster-Zuber [21] are also 
presented for comparison, The data for full surface nucleate boiling 
agrees well, Note that the surface temperature remains approximately 
constant until full nucleation is established. The reproducibility of 
this surface is shown in Figure 13. Arbitrarily, mirror run #2 was 
chosen for all subsequent comparisons. 

_— ye Effect of Roughness 

Figure 14 shows, comparatively, the effects of surface roughness of 
the copper test sample. Roughness affected not only incipient boiling 
flux, but wall superheat temperature as well (see Table I). As expected 
increased roughness provided more nucleation sites at lower wall super- 
heat, 

The effect of the Teflon coat on the roughened surface at low heat 
fluxes was that of partial insulation, Since nitrogen wets the T2flon, 
no non-wetting interfaces occurred such as observed by Young and liumel 
[26 | with a water-—Teflon-metal combination, which produced heat transfer ; 
coefficients ten times the plain roughened surface:(Fig, 19). However, the 
Teflon provided areas of low free surface energy (whereas the copper had 
a high free surface energy {27} ), altering the surface tension forces 
and allowing a somewhat higher heat flux for a given AT, 

A plausible explanation of the Teflon curve crossing the rough 
surface curve is that, originally, a greater surface temperature was 
required to activate a site located on the Teflon (bubble nucleation 


point seen visually) due to the insulating effects of Teflon (similar to 
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grease coating discussed in the following section), but that, once 
energized, greater heat transfer could take place due to the greater 
number of activation points, The bubbles formed were smaller in diameter, 
but larger numbers of them were seen per surface area, giving an overall 
higher heat flux. 
Effect of Contaminants 

The effect of surface contaminants (Figure 15) is that of insulation 
requiring higher surface temperature for a given flux, This effect was 
expected for the oxidized surface [ 5] as well as the greased surface[ 1] 3 
but, because of an insufficiently thick layer of oxide, the boiling 
curve for this surface coincided with the results of the mirror finish- 
un-oxidized. It was not possible tore-investigate the effects of surface 
oxide thickness for this study, 

The surface coated with Apiczon "N" grease followed expectations 
with increasing heat flux. Heat flux change was approximately proportional 
to QT up to 5 watts/cm*, indicating a conduction layer, At full surface 
nucleation (5 watts/em*), the slope assuwned the characteristic nucleate 
boiling form. Boiling, however, causes a physical change inthe grease 
surface, discussed in greater detail in the following section. 

Hysteresis Effects 

The differences in observed superheat when the heat flux is increased, 
reduced, and again increased, is called nysteresis. This efteae is 
vividly demonstrated in Figure 16, showing the Teflon and creased surfaces, 
It was also observed to a much smaller extent (not shown) in the mirror 
copper surface, The explanation is that for the original flux increase, 
a few sniall nucleate patches became activated. At increasingly higher 
heat fluxes these patches spread, activating nearby sites. Upon decreasing 


heater power, thcse active sites were observed to thin out evenly, but 


30 


many remained active at fluxes which previously had caused only convection 
heat transfer or partial boiling. Since the time for a data taking run 
was usually an hour, changes in surface entrapped vapor were expected. 
As the surface boiled, more noncondensable vapor was boiled off. The 
final effect was that more energy and consequently a higher AT was 
required to reactivate a surface as the population of any one sized 
fertile site (one containing vapor) had been reduced. The Teflon coated 
sample followed expectations. Observation of this surface showed that 
for a small heat flux (.4 watts/cm2), a small area of sites was activated, 
With an increase in power, more of the surface was covered with nucleation 
sites, The number of active sites increased rather slowly. It was not 
until a heat flux of 12 watts /cm2 was established that the whole surface 
became covered. On decreasing flux, all sites remained active, however, 
and at 0.4 watts/cm’, the bubble columns were very fine, hair-like streams, 
indicating many active sites. At this condition, the lowest AT (.3°K) 
was recorded for nucleate boiling with any surface considered in this 
investigation. The reactivation wall temperature was higher, as expected, 
The grease coated surface showed a wider hysteresis curve because 
of the aforementioned altered surface condition. This change was noted 
on power reduction. At each point where nucleation had occurred,a small 
indentation appeared in the grease, These indentations were distributed 
evenly over the surface. The plausible explanation is that the surface 
became superheated enough to transfer energy in such great quantity to 
the entrapped vapor on the polished surface thereby activating nucleation 
sites. Bubbles formed from these sites with enough force to deform and 
remove the grease, (which acts as a brittle solid at liquid nitrogen 
temperature), Differing from the Teflon surface, however, it was noted 


that the reactivation energy was lower than activation energy because 
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the surface resistance to nucleation had been decreased due to the "pits't, 
It was intended to record this phenomenon photographically in order to 
provide data on the number of active sites per unit area. Upon dis- 
mantling and warming the surface to ambient temperature to get a close- 
up photograph, however, the grease was observed to flow partially thus 
diminishing the effect. 

Effect of Materials 

The effect of material is shown in Figure 17, The results show 
better heat transfer with copper than with nickel. This result agrees 
with Berenson's | 5] results for pentane, The data for platinum [13] 
is also shown, The thermal conductivity of platinum is very similar to 
that of nickel and, as can be seen when the surface is nucleating fully, 
the curves are very similar, 

Due to the lower thermal diffusity (<= wo ) of nickel compared 

= Ww 
to copper at LN» temperatures, the bubble generation rate is lower, 
resulting in poorer heat transfer. A plot of < versus wall superheat 
(Figure 18) clearly shows that as the thermal diffusivity of the boiler 
surface increases, the wall superheat decreases (and hence higher heat 
transfer coefficients). 

The nickel surface was unusual in that it was the only surface 
investigated which exhibited a temperature drop for a heat flux 
increase (Figure 17). As yet, no explanation exists for this negative 
slope phenomenon [24] « There may exist an unstable condition, as 
observed by Marto [14] with sodium pool boiling. The basis for this 
conclusion is that the thermocouple closest to the surface was observed 
to flucturate = 34V (0.3 K overall) at 5.8 watts/cm, finally settling 
to + 1,V at 13.3 watts/cm2., Virtually no hysteresis occurred on power 


a 


decrease, and this same unstable thermocouple fluctuation was observed 
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in the vicinity of 10.4 watts/cm@ heat flux, 

A general comparison covering the spectrum of materials and fluids 
used in heat transfer equipment is shown in Figure 19. Of interest is 
that, in general, the observed heat transfer coefficient is very nearly 
5 x 103 BTU/(hr ft? °F) for all cases at the highest observed tempera- 


tures (fluxes). 


33 


SECTION 5 
Conclusions with Recommendations 

The results discussed lead to the following conclusions: 

1. Surface conditions significantly affect the incipient boiling, 
the nucleate boiling wall superheat and slope of the boiling curve of 
liquid nitrogen, 

2. Boiler history (i.e. whether previous boiling has occurred on 
the surface and whether the heat flux is being increased or decreased) 
affects the nucleate boiling wall superheat, (T,-T,). 

3. Boiler surface materials affect the boiling curve of liquid 
nitrogen. 

Therefore, correlating equations such as that of Rohsneow and 
Forster-Zuber, should be used only when all surface conditions are specified 
and account has been made for each variation from a reference condition. 

The results indicate that further studies be continued in order to: 

1. Completely eliminate all "rogue"! nucleation sites on the boiler 
test surface. 

2. Determine the effect of surface conditions on bubble frequency 
using a high speed motion picture camera, 

3. Continue studies with heavier layers of oxide and grease to 
further bear cut the effects of contaminants and add further credence 
to application of the entrapped vapor phase nucleation site theory at 
cryogenic temperatures 3 | ° 

4, Study the effect of wettability and non-wettability in conjunction 
with small sized or artificial cavities. 

5. Widen the realm of cryogenic fluid heat transfer studies by 
using liquid helium as a working fluid and investigate the effect of 


surface conditions. 
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APPENDIX A 
Calibration 

Because a high degree of accuracy was desired, direct instrument 
calibration against secondary standards was chosen as the best means to 
attain the results. The amplifier-voltmeter combination was set up in 
an electonically noise-free room, using regulated power in order to 
reduce electronic disturbances, A standard cell of 1.01944 volts was 
used in combination with a voltage divider to generate small signals 
(0-600 aN ). The amplifier was then adjusted so that the output displayed 
on the digital voltmeter agreed with the value determined by a Rubicon 
potentiometer known to be accurate to within SyAVe Periodically, for 
the next three days, amplifier drift was checked. It was determined to 
be less than Jan No change was observed using unregulated power. 
A month Water, calibration was re-checked against a Leeds and Northrup 
K-3 Guarded Precision Potentiometer, The amplifier-digital voltmeter 
combination agreed to + ea over the Oa: The accuracy of the 
K-3 potentiometer is 2M ° 

Using techniques outlined by the National Bureau of Standards [17] ‘ 
the thermocouples were checked for homogeneity after installation. Due 
to the limited range over which the thermocouples were to be used, cali- 
bration was accomplished by using two known boiling points covering the 
temperature spectrum of expected data points. In each case, the reference 
used was liquid nitrogen, Aviator’s liquified oxygen (99.95% purity) 
was used as the prime calibrating point. Data was recorded for the LN,- 
LN» combination prior to LNo- LOX. The small zero readings were subtracted 
from the liquid oxygen reading for each T. C. The known boiling points 


were corrected (method in Appendix C) to standard pressure (760. mm Hg). 
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The LOX boiling point was extremely stable, Generated emfs remained 
within Pe over a half hour period and were reproducible on a 
succeeding day. 

This data then became the input to a computer program developed by 
the NBS Cryogenic Laboratories at Boulder, Colorado [16] » The progran, 
described fully in Appendix C, compared the resultant spot calibration 
with the NBS calibration table, calculated a correction factor, then 
generated a working table for each thermocouple, An option was built 
into the program so that the printed output tabulation of temperature, 
emf and slope was based on a liquid nitrogen reference. An example of 
the output is presented in Appendix C also. 

In order to check the accuracy, the boiling point of liquid methane 
was used as a check. Ultra pure methane (99.95%) was liquified in the 
laboratories. For an as yet unexplained reason, the temperature of the 
LCH4 decreased with time. It was assumed that the LCH4 was undergoing 
a contamination with the atmosphere. Notwithstanding, an emf-time graph 
was drawn for each thermocouple tested over a period of six minutes, and 
the subsequent straight line was extrapolated to a zgero time value 
(completion of liquification). The time between liquification and first 
datum point was 30 seconds. This temperature decay with time was obser- 
ved each time the methane was liquified. The extrapolated values were 
used to generate working tables for each thermocouple, The methane 
tables were compared with the oxygen-based tables and agreed within 
364 pV (.16°K) at liquid methane temperature (111.6°K). ss 

A summary of both generated working tables at several arbitrary 


temperatures appears in Table A-l. 
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TABLE A-] 
WORKING THERMOCOUPLE TABLE SUMMARY 


Difference 
_LCH, ,V{ _LOX-LCH, V 





Diff. 
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APPENDIX B 
Generation of Thermocouple Tables 
Working thermocouple tables were generated using guidelines set 
forth in NBS 8750 [16] and program FACTOR with standard copper-constantan 
(TP-TN) thermocouple table decks. This program with tables was available 
at cost from NBS Cryogenic Laboratories, Boulder Colorado. It was 
written in FORTRAN II, but modified by the author to conform to FORTRAN 
60 for submission to U.S. Naval Postgraduate School Computer Facility 
using CDC 1604. This modified version of the program is presented along 
with the standard table and input data at the end of this appendix. In 
addition, selected cards were pulled from the program deck and collected 
together in a manner suitable for checking spacing of card information 
and commands. Finally, an example of a generated working table has been 
included. 

The computer compares the observed temperature change in degrees 

a oe with the standard table values and calculates a factor from this 
ratio Cnv sagen cdg standard). A table is then printed from 0 to 
300°K for any selected reference temperature in this range. 

Input data to the program appears on two or more cards depending 

on how many different tables are being processed sequentially on the one 
-program submission. Card A contains calibration information, the details 
of which are as follows: 

Column 1-63 Temperature scales; either DEG. K or DEG. C. Units must 
be consistent throughout this card. The space between 
the period and K or € is required. 

Column 11-203 High Temperature of spot calibration which must include 
a decimal point. 


Column 21-302 Low temperature of spot calibration with decimal point, 


Al 


Column 31-403 Voltage difference observed from spot calibration 

in microvolts with decimal point. 

Column 41-503 Reference temperature desired for tables with decimal 

point, 

Column 58-60: Punch YES if punched cards desired with table-- 

otherwise punch NO in columns 58-59. 

Card B is used as an identification card, the details of which are 
as follows: 

Column 1-153: Company or laboratory identification 

Column 21-263 Six characteristics available for lot number, or 

spool identification, 

Column 31-423 User’s name, 

Column 51-643 Date of test. 

Card C is used as a repeat card and must preceed subsequent to A and B 
cards for multiple table output. It should be exactly as follows in all 
casess 

Column 1-213 Punch REPEAT LAST T. C. DECK. 

Column 503 #£Punch 0 (Zero). 

Column 55-603: Punch DEG, K with space between period and K, Use 

DEG. K regardless of units wanted. 

Column 68-703 Punch YES. 

Assembly of the deck for submission to the computer from top to 
bottom wass 

a. Computer monitor control cards 

b. FACTOR program deck 

c. Standard data deck 


d. Card A 
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e. Card B 

f. Repeat card 

ge Card A 

h. Card B Multiple table output 

i. Repeat card 

lee “G@BCo sc 0c 00 

Computer run time for programing 18 sets of data was two and one-half 


minutes, 
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APPENDIX C 
Sample Calculation with Error Analysis 

The following are assumptions used in reducing data: 

a. Temperature tables in microvolts are accurate to within = 
(see Appendix A) from 77.3 K to 90.1°K and within ¢ 2 pv above 90.1°K, 

b. Liquid nitrogen used had a minimum purity of 99.95% and the 
same quality used for calibration. 

ce Thermocouple hot junction positions are known to t .O1 inches, 

d. Thermal conductivities of electrical tough pitch copper agree 
with tabulated values within 5%. 

e. Barometric pressure was accurate to + 0.1] mm Hg, 

f. Measured depths of nitrogen are to within 1 cm. The reference 
thermocouple was always in 10 cm of nitrogen, 

Arbitrarily choosing point 8 of the oxidized copper surface with a 


barometric pressure of 766.8 mm Hg at 21.2°C. 


TABLE C~1l 
ZERO CORRECTED AND REDUCED THERMOCOUPLE DATA 
T. Ce Location from Boiling Corrected Temperature 


Surface (inches) for Barom, and Reference 
Immer, Depth °K 


1 84.42 = .12 
2 85.90 * .12 
3 Siesus oll 
4 88.90 = .11 





The saturation temperature correction for barometer and reference 
immersion depth was similar to the calculation of the saturation temperature 


correction used for working fluid shown below, 
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T Saturation 
Depth of LN» at beginning of run was 42 cm and at the end was 32 cn, 
therefore, for error purposes, the average was assumed to be 3725 cm, 


Using the standard pressure depth equation and changing units: 


eS ru (32.Alb) (3225¢m \ f+ > ) (2222 a 
e 2 =i--:*=<“(<i‘(C«];t;CS*~*;: Ts a ee —- =< = 
. ft3 /N\2 54emAn \I728in4 \ 14.7Ib/iné 


Phead = 21.96 + 2.97 mm Hg if 


APh = Ph — Po = 760.0 - 766.80 +.1 (C-2) 
APp 


WW 


6.80 +.1 mm Hg 
Fo wae taken by four way interpolation from Ref, EY for barometric 
pressure and temperature, 
P, = -2,.3mm Hg 
APeotal = Phead + APL +P, (C-3) 
= 21.96 + 6.8 - 2.3 


26.46 + 2.97 


OProtal 


Making use of the perfect gas relationship, the Clausius-—Clapeyron 


equation takes the form: 





at, «Rls SFotal (c-4) 
heg Pe 
oy Ib (139 °R) @ (Phot mm ra) (24 
Ib 85.7 btu 760.mm hg/\9°R 
= .1167APtota] Kk 
Substituting from aboves 
AT, = (,1167) (26.46+ 2.97) 
AT. = .309+ .035 K 
T. = T, (tabulated [27] )+ aT, (c-5) 


= 77.347 + (.309 + .035) 


-~ 77.66 + .04°K 
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Heat Flux 

By plotting T. C. temperature readings versus position from the 
boiling surface with error limits included, and connecting the extreme 
points of the data, minimum and maximum sloped lines (4T/X) for these 
points can be determined, 

Using the Fourier conduction equation and knowing the thermal conduc- 
tivity k, based on average temperature (°K) between T. C. position 1 and 


4 (Figure 10), 











Keating Karam (=<) a In. oa 


A ax NemeK/ \7o3m/ \a5acm 
For the sample data listed, Q/A maximum = 12.42 watts/cm@ and Q/A minimum = 
11.20 watts/cm2 for an average of Q/A = 12.31 ¢ 1.1) watts/cm@, 
T Wall 
Altering equation (C-6) gives: 


T, - Th = 0/A = (C-7) 


where subscript 1] indicates thermocouple 1 position. 


Substitution of the average Q/A into (C~7) gives: 


T) - Ty = (12.3iw)(.06 inches) (2.54 cm) (cm°K ) 
Ge in) 5.03" 
T) - Ty = -373°K 


Using the value of T) = 84.42 from table B-l 
Ty = 84.42 - .373 = 84.05 °K 


The most probable error for the wall temperature is; 
Ut. (amp 4 (Ty = 79) aa ae ced (C-8) 
Q/A K Xy 


fractional error in heat flux = ,.0891 


eB 

ie 

> 
ij 


fractional error in thermal conductivity = .050 
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fractional error in thermocouple ] position = .166 


4T, = + .1K° from the temperature reduction (Table B-1]) 
Substituting in equation (C-8) 
AT, = + .14°K 
Twall - Tsaturation 
Simple subtraction gives the following for a value of 
AT =T, -T, = (84.05 t 0.14) - (77.66 + 0.04) = 6.39 + .15.K 


. 
i 4, 
oe ie 
Gomes 
Aged 
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FIG, 16 HYSTERESIS EFFECTS WITH TEFLON AND APIEZON ''N'' GREASE ON COPPER SURFACE 
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